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Abstract 

We elaborate on the computation of the pressure of thermal quantum electrody- 
namics, with massless electrons, to the fifth (e^) order. The calculation is performed 
within the Feynman gauge and the imaginary-time formalism is employed. For the 
calculation, the method of Sudakov decomposition is used to evaluate some ultravio- 
let finite integrals which have a collinear singularity. For the contribution, we give 
an alternative derivation and extend the discussion to massive electrons and nonzero 
chemical potential. Comments are made on expected similarities and differences for 
prospective three- loop calculations in QCD. 
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1 Introduction 



Recently we reported on the order and contributions to the pressure of 
a massless QED plasma at temperature T and zero chemical potential fie- Here we 
would like to fill in the discussion, particularly for the more difficult calculation, 
as some of the techniques might be of wider interest. In addition, an alternative 
derivation of the contribution will be given. 

For motivation we note that the equation-of-state of a relativistic plasma is of 
relevance in astrophysics 0, ^, ^ ^ . It was obtained at nonzero T and fie by Akhiezer 
and Peletminskii |Q to the third order (e'^) and later extended to order at T = 0, 
but nonzero fie, in both QED and quantum chromodynamics (QCD)[^, |]. Thus the 
calculation of the T 7^ 0, /ig = contribution at order fills a gap in our knowledge. 
It is also a step towards same order calculations for QCD which is believed to exist in 
a perturbative quark-gluon phase at high temperature. Currently the pressure of the 
high-temperature phase of QCD is known to order g'^lng, excluding the normalisation 
of the logarithm |05 H- 

The calculations in this paper are performed within the framework of the imaginary- 



time formalism (see |]T0|, [TT| and references therein) whereby the energies take on 
discrete Matsubara values, go = iniiT, n being an even (odd) integer for bosons 
(fermions). The usual zero-temperature ultraviolet (UV) singularities are regularised 
by dimensional continuation [|12| . For simplicity we restrict ourselves to the Feynman 



gauge but renormalization via minimal subtraction ensures that the coupling 
constant is gauge- fixing independent |]T^|, and hence so will then be our final answer 
for the pressure |Q. 

At intermediate stages of the calculation we will encounter various types of infrared 
(IR) singularities. The first kind is due to many-body effects and gives rise to power- 
like singularities in some diagrams. As is well known, these are removed when the 
static electric propagator is dressed to take into account the screening of electric 
fields in a plasma [^, ^ |TT|. As a result of this resummation the expansion for the 
pressure is in powers of y/^ rather than e^, the famous "plasmon" of Gell-Mann 
and Brueckner |]l5l illustrating this at the lowest order. In this paper we will obtain 
also the order plasmon contribution. 

We remind the reader that the identification of bubble (i.e. no external legs) 
diagrams which require the use of screened propagators so as to produce a consistent 
perturbative expansion is most easily done in the imaginary-time formalism where 
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only the zero mode of the photon propagator lacks an infrared cutoff of order T. 
More generally one can perform a consistent resummation in imaginary-time for any 
Greens function whose external legs are bosonic and static (at zero energy). For 
Greens functions which have external fermionic lines, or nonstatic external bosonic 



lines, one must first analytically continue |16, 11 1 to real-time to obtain the physical 



Greens function before the resummation can be discussed. In the latter case it is in 
general necessary to use nontrivial propagators and vertices to restore the perturbative 



expansion, as discussed by Braaten and Pisarski |jT7|. 

Though individual bubble diagrams in imaginary time may be IR finite (after 
using dressed propagators if need be), once the frequency sums are performed the 
diagram in general splits into several pieces (integrals) each of which individually may 
contain mass-shell and/or coUinear singularities. Since it is convenient to evaluate 
the different integrals separately, a regularisation has to be used for these "spurious" 



IR divergences. We will use dimensional regularisation (DR) for this also [T^, ^ 
^. Actually, in this paper we will make no attempt to distinguish between IR and 
UV divergences at each stage but will only verify that sums of diagrams which a 
priori should be finite, are indeed so. Such finite sums form gauge-invariant subsets, 
contributing to the pressure amounts which are conveniently labelled by different 
powers of A^, the number of fermion flavours. 

We mention here also two other difficulties, related to IR divegences, in evaluating 
some integrals. As in the case of the 3-loop pressure in 0^ theory |l^, we find some 



integrals containing singularities along the path of integration which are intrepreted 
in the principal value sense. In some cases we will transform, by a change of variables 
in the integral split at the point of singularity, these principal value singularities into 
infrared singularities and then evaluate them by dimensional regularisation. Also 
some difficult 3-loop integrals which are UV finite because of statistical distribution 
factors, but contain a collinear singularity, are in this paper handled by using Sudakov 
variables in dimensional regularisation (see pT|). 

Our notation is as follows : the wave-vector, = (go, q) is contracted with 
a Minkowski metric, = Qq — (f, and the measure of loop integrals is denoted 
compactly by 



/ 



[dQ]^T J 



qo,even 



(27r)^-i 
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/ {dQ} ^tj: I 



go, odd ' 



d 



D-l, 



D-l 



The fermions are kept as four- component objects, Tr^'y^'j^) = ig^^,, and the gauge 
propagator is D^y{K) = g^y/K"^. 

In the next section we review the calculation of the pressure in massless QED 
to third (e'^) order. In Sect. 3 the diagrams and integrals which contribute at fourth 
order are discussed. Some frequency-sums are evaluated in Sect. 4 while in Sect. 5 we 
explain our use of Sudakov variables to evaluate two difficult integrals. The final 
pressure to order is summarised in Sect. 6 and the renormalisation group briefiy 
discussed. In Sect. 7 we rederive the contribution to the pressure. The results of 
this paper are summarised and discussed in Sect. 8 while the appendices contain some 
useful identities and technical derivations. 

2 Lower Orders 



Before discussing the 3-loop calculation, let us briefiy review the lower order results 
for the pressure of QED with massless Dirac fermions. The ideal gas pressure Pq 
due to electrons, positrons and photons is given by, 

T 



Pn 



V 



In 



Det+{d'g^^) ^ Det+{d') Detail ^ 



7 



{D-2) + -{AN) 



(2.1) 



TT 

45 



i + \n 



In the ± subscripts refer respectively to periodic and antiperiodic boundary 

conditions, and the second determinant in is the ghost contribution which is required 
for proper counting of physical degrees of freedom in the ideal gas pressure |T^ . The 
first correction P2 is given by the two- loop diagram of Fig.l, 



I 



K'^{K-py 
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= (D-2)e2iV/i(26i-/i)/x^-^T2^-^ (2.2) 
288 

The last line above follows from the previous one as D — > 4. In the rest of this paper 
we will often take this limit, where consistently possible, without comment. In ( |2.2| ), 
e is the renormalised coupling, /x the mass-scale of dimensional regularisation and the 
integrals bi and /i are defined through 



fn — 



[dQ] 

(g2)n 
{dQ} 

{Q 



2\n 



We have scaled out the temperature in the above integrals so that T = 1 there. The 
simplest way to evaluate 6„ is to first perform the momentum integrals and then the 



frequency sums [Q. This gives 



(27r)2'^ T{n) ''^ ^ ' 



Then is easily obtained by using a scaling argument : Consider the sum /„ + &„ 
and scale the momenta in the integrals by a factor of 2. Thus fn + bn = (2^"^^"^)6ri, 
and so 



U = (2^"+i-^-l)6„. 

The next contribution to the pressure is of order (e^)^/^ and reflects Debye screen- 
ing [0, 1^, ^, 1^, |Tl|. To lowest order the static electric polarization tensor noo(0, q 
0) = m? + 0{q^) with m = eTN^^"^ / ^/S. Insertions of noo(0, 0) along the photon line 
in Fig.l create IR divergences which sum to 



2 J {27rr^, p \q 
T f d^q 



1 1 
m , m 



ln(l + -^) - (2.4) 



2 J (27r)3 V 



4 



- W(f) • f^-^) 

Though ( p.4[ ) is finite, we have evaluated it using dimensional continuation 3 
D — 1. Then the second term in ( |2.4| ) vanishes and the first term gives the result ( |2.5| ) 
as D — > 4. The point of this excursion is to verify, in this example, that scaleless 
integrals like those above may indeed be consistently dropped in DR. This fact will 
be useful later. 

It is of interest to note that one may also calculate the plasmon contribution in 
real-time but then the analysis is much more intricate. In a real-time analysis, the 
transverse photons are relevant at intermediate stages but their net contribution to 
the plasmon term vanishes [^4|, and one recovers the result of the imaginary- 
time analysis where it is clear from the outset (by power counting) that only the 
longitudinal photons contribute. In the language of Braaten and Pisarski [O , for the 



calculation of the pressure of QED in imaginary time (recall the discussion in Sect. 
1), the only soft line is the static photon line. Thus resummation (which gives the 
e^, etc. terms) involves dressing this soft line with the relevant static "hard thermal 
loop" which is just the electric mass noo(0,0) [jl^- For more discussion comparing 
the real and imaginary-time approaches to the calculation of the pressure see [53, 



in particular the conclusion of the second reference of |24 



3 Three Loop Diagrammatics 

The order e^N diagrams are shown in Fig. 2 and 3. For massless fermions, the Ward 
identity Zi = Z2 implies the mutual cancellation of the counterterm diagrams (Fig. 3), 
so the sum Gi + G2 (Fig. 2) is UV finite. After performing the spinor traces and 
some algebra. 



Gi = T'""-' /-'^ 2 {D - 2f (Uh - h,f + h- 2/5) , (3.1) 



G2 = T'""-' /i«-2D 2 (2_D) ((861 - 16/i)/i + 2{D - 4)/2 + S{D - 3)h + ^h) 

(3.2) 



5 



As before, we have scaled all the momenta by 1/T so that the integrals are dimen- 
sionless (i.e. T = 1 there). The integrals 6„ and /„ were defined in the last section, 
while the rest are 



, {dKdR} 

{dK}[dQdP] 
K^Q^P\K + Q + Py ' 

{dKdRdS} 

{dK}[dPdQ] p2 
K^Q^{K + Q)\K + P)\K + P + Qy' 

{dK}[dQdP] (P-Q) 

K^p^Q^{K + Qy{K + py ■ 

Some simplification is possible. Firstly, within DR, Ji = as shown by Arnold 



and Espinoza |23] using scaling arguments. We will deduce the same result by direct 



evaluation in the next section. Using scaling arguments as in ||23[ one may also show 



i (2"-^^ - l) i/i - Ih (3.3) 



where 

[dQdPdK] 



Hi 



K^Q^P^{K + Q + Pf 

is the integral analysed by Frenkel, Saa and Taylor [|1^ in their 3-loop pressure cal- 
culation in hot 0^ theory. Next, some algebraic rewriting gives for J4, 

{dK dRdS} {K + Sy 



Km^S^K + Ry{K + R + S) 



2 /i /i - Ih 



-I" 



2 



The last line is valid because /i is finite near D = 4 but Ji = as mentioned above. 
Relabelling J3 — H2 and J5 ^ H^, we obtain 



6 



Gi + G2 



{D-2) 



e4 ]V ^&-2D 



(2(1- 2("-3^)) Hi + (20 - 3L>) H2 



6 

(D-2)2(2/f3 - /2(/i-M')- 



(3.4) 



Consider now the order e^N"^ diagrams. The photon wave-function renormahza- 
tion required for (Fig. 4a) is provided by Xi (Fig. 4b). Diagram G^ also has an 
IR singularity which contributes to the plasmon as the first term of the series in 
eq. (|2.3|) . In principle this term should be subtracted from G^ to avoid overcounting 
but since, as discussed earlier, it vanishes in DR, double-counting is automatically 
avoided. We have 



G, 



+ , (3.5) 



and 



T 



A-D 



Xi = -(Z3 - 1) iV (D - 2) T^^-s h (26i - /i) ) , (3.6) 



where 



N 



+ 0{e% 



67v^{D - 4) 

In ( p.5|) we have dropped a vanishing contribution proportional to Ji . The new integral 



m 



IS 



[dQ]{dKdR} {K-Rf 
Q^Km^{Q + Kf{Q + Ry 



In summary, our task for the calculation of the order contribution to the presu- 
ure as given by eqns. ( |3.4| , P75| , [H^ ) has been reduced to the evaluation of the following 
four integrals : 



Hi 



Ho 



[dQdPdK] 



{dKdRdS} 
K^R^S^{K + R+Sy ' 



(3.7) 
(3.8) 



7 



rr ^ f {dK}[dQdP]{P.Q) 

^ J K^p^Q^K + Qy{K + py' ^ ■ ' 

^ r [dQ]{dKdR}{K-Rf 
" J Q^Km^Q + KfiQ + Rf ■ ^ ■ ^ 

In the next section we will do the frequency sum in Hi and re-obtain the expression 
of . The analysis of H2 is completely analogous to Hi . The new integrals we have 
to analyse are i^s and H^. These are the main concern of the next two sections. 

4 The Frequency Sums 



There are several variations in the literature |]T^, ^ for performing the frequency 
sums. Here we sketch one way (cf. [p6[] ). 

•h 

Consider the two-loop integral 



The first step is standard and involves writing the sum as a contour integral, 



1 _ r d^-^ K 1 r -Nk, 
where 

1 



fco 



and the anticlockwise contour circles the poles of N^^ only. Next the contour is 
deformed to enclose only the poles of the propagator. Upon using the identity 

1 r°° 

— {f{x)-f{-x))= / dxo6{xl-xy{xo)f{xo), (4.2) 

ZX J —00 

and noting that go is an odd multiple of i, the result can be written compactly as 
d^ K , , r 1 



where 
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The above procedure is repeated on the go sum but note that ko is now real and 
continuous, 

f d^ Kd^Q ejko) 6{K') efa) 6{Q^) 
^1 - 7 ^27r)2^-2 (ZTqP ^ " ~ ^ ~ '^""^o^ ■ 

Using the identities in the appendix and symmetrizing, this becomes 

Now use the identities in the appendix again to write the expression above in terms 
of rifc and Nq where k = \ko\ and q = \qo\- Then terms which are independent of 
statistical factors contribute zero in DR, terms with only one statistical factor also 
vanish by Lorentz covariance and DR, while terms with two statistical factors vanish 
in DR after the energy integrals are done. Thus Ji = in agreement with . 



•Hi and H2 

When applied to Hi, the frequency sum algorithm discussed above yields 
[dK dQ dP] 



Hi 



K^Q^P^{K +Q +Py 

d^K d^Q dDPe{ko)6{K^) ejqp) 6{Q') ejpo) 6{P^) 
(27r)3(^-i) (K + p + qy 

d^K d^Q d^P6{K^) 6{Q^) 6{P^){4:nkngnp + GukUp) 



(27r)3{^-i) iK + P + QY ^^'^^ 

where k = |/co| and we have skipped a few steps which make use of the identities in 
the appendix. We have also dropped terms which vanish in DR. To write (|4.3| ) in the 
form of Ref . UTOl , use 



r r°° 6(K'^)6(P'^)6(Q'^) 1 /" f°° 

jd^-\j_Jkodq,dpo ^K^Q^py - —Jd^'^qjjkodqodpo 

5{K^)5{P^) 



(g2_^o+)((K + g + p)2-2 0+) ' 

(4.4) 



to get 



-6 fd^Kd^Qd^P 6{K^)6{P^)nkn^ 
Hi = 7r~ 



2m J (27r)3(^-i) {Q"^ - i()+){{K + Q + Pf - 



r d'^K d'^Q d'^P 5{K')5{Q')5{P') rip 
J (27r)3(^-i) {K + Q + Py ■ ^ ' ' 
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This is the expression of |]T9[. Keeping only the nonvanishing terms as D — 4 gives 

m 



Hi = , \ ( 5.6658 - — -I . (4.6) 



We leave it as an exercise for the interested reader to show that ( |3.8|) is obtained 
from (|4.5|) by replacing nj with — A^j, and then 

H2 = , \ I 2.6045 - — ^ — - I . (4.7) 

3(2%2) \^ (^ - 4)7 ^ ^ 



For we obtain 



^3 = / '^''^^^'^^Xf^ 6(go) 5(Q^)e(A:o) 5(i^^)e(po) 5(P^) {7^}, (4.8) 
where 



J-7?l = I (TxliT^i^-Q) {Tk\K^K-P) 



(ir + Q)2(ir + P)2 {^K ~ Q)\K - Q \ Pf - P)\K - P ^ QY' 

Tk = Nk,{npo[ngoP-Q-N^goP-iQ + K)] 

-N_,,M -(P + K)- + K)-{P + K)]) , 

and {Tk\K K — Q) means replace by — in the expression for Tk- After 
the usual simplification, 

H3 = Ji + Ki + L,. (4.9) 

The piece Li contains integrals that can be performed exactly using relations 
(|AT7|)- (|A:28|) of Appendix A. 



Li = \^^y\D-i) 2(^-^) 5(1/2, D/2 - 1) B{D/2, D/2 - 2) 
xM2{D-5) [Mi{D - 3) + M2{D - 3)f 
, ..(D).(D_1) ^^^^2 - 1, D-3) 2-- AT,, 



(27r)3(^-i) 
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where 



+2 



(T,7=±l 

. , , - 10\ , , /SD - 10 
A^i ^ + A^2 



E C'i(a,7). 



The functions cj, Aii and are defined in the appendix and B{a, b) is the usual 
beta-function. Expanding about = 4, 



-1 /27r2 + 77r2ln2 + 67r2C'(0) + 18C'(2) 



3(2%4 



P-4) 



+ 35.478 



(4.10) 



In ( [4. 91 ) i^i is an integral similar to that appearing in Hi and H2, 

d^K d^Q d^P 6{P^)6{Q^)6{K'^)NpNgnk 



(2vr) 



{K + p + qy 



1.1439 



2(27r)e 
Finally, Ji is defined by 



(4.11) 



Ji = J "^""l^jirclf ^ '^+(-^')'^+(^')'^+(-P')^^o^.o(iVpo + n,,) E (-^)'5i(a,7), 

(4.12) 



with 



'5i(a,7) 



P-Q 



1 



• Q • Q + P ■ (ai^ + 7Q) ' 
and 5+((5^) = 5(<5^)6'(g°). The evaluation of Ji is described in Sect. 5. 
Collecting the pieces, 

Hs = Ji + Ki + Li 



(4.13) 



1 0.4112335165\ 
1.095 — 



277r4 



(D-4) ; 



(4.14) 



The only extra point here is the doubled propagator 1 / (Q^)^. Thus the contribution to 



11 



the contour integrals coming from this propagator must be extracted correctly using 
the appropriate formula of complex analysis. Other than this technical point, the rest 
of the derivation proceeds as discussed for but the expressions are lenghtier. So 
we will only state the results. 
We obtain 

H^ = J2 + K2 + L2, (4.15) 
where J2 is an integral similar to Ji and is discussed in the next section. 

d^K d^Q d^P 5{K^)5{Q^)5{P^) Ng Np 



(27r)9 (K + p + qy 

0.4417 



and 



4(27r)6 ' 

/ 177r2 + 127r2 In 2 + 247r2C'(0) + 72C'(2) 



33 210 7r^\ {D- 4) 

Finally 

Ha = J2 + K2 + L2 



0.5098 



2V (D-i) ) 



(4.16) 



(4.17) 



1 /,.,o,_ 2.254840072\ 



5 Sudakov Variables 

In this section we describe the evaluation of the integrals Ji and J2 appearing in the 
last section. 
•Ji 

The integral Ji was defined as ( [4.12| ), 

X / d^'K 5+{K^) ^—^ ^ ^ . (5.1) 

J +^ ' K-QK-Q + P-{aK + ^Q) ^ ' 
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This integral is quite complicated and we have not succeeded in evaluating it in closed 
form. However for our purposes we only require the order {D — 4)~^ and 0{D — 4)° 
terms from (|0| ) in the hmit D ^ A. Although ( ^.11) is UV finite, a l/(i^ — 4) pole 
is expected because of the collinear singularity as K.Q 0. 



Let us concentrate on the i^'-subintegral appearing in (|0 



For the rest of this subsection we will mostly discuss this subintegral with the implicit 
understanding that it occurs inside so that the constraints = = and 
Po = P, Qo = Q hold. It is sufficient to extract only the 0{D — 4)'^ and 0{D — 4)^ 
pieces from (|5.2| ) because the remaining (P, Q) integrals in (|5.1| ) are UV and IR finite 
( so that we will only need their 0{D — 4)° and 0{D — 4)^ pieces to get the full 
0{D — 4)"^ and 0{D — 4)° terms for Ji when our result for ( |5.2| ) is substituted back 
into dni) ). 

Since the integral ( p.2| ) involves three independent scalar products in D dimen- 
sions, it poses a formidable problem. We begin our task by decomposing the loop 
momentum in a Sudakov base constructed on P and Q in such a way that the 
D dimensional angular variables will disappear from the denominator in ( |5.2| ). Thus 
we write 

K = aP + xQ + Ks_ (5.3) 

with K± denoting a transversal vector (P ■ K± = Q ■ = 0) oi D — 2 nonzero 
components. Note that the delta and step-functions in ( 4.12| ) imply Kj_ < 0. Since 
the Jacobian in going from the i^'— basis to the (a, x, K±) basis is s = P.Q, we obtain 
for ( |5.2| ) (upon using ( |5.3| ) to eliminate P.K and Q.K ), 

~ da r I d^-'K, H2sxa + Ki) N, 9(k,) 

' J a a + (cTx + 7) ^ ^ 



00 J —00 



- d. I ,^-'^k' (5.5) 



where now 



ko=^-l^P + xq + iK^)o , (5.6) 
2sx 

which follows from (|5.3|) and the fact that po = p and Qo = q. The original collinear 
singularity of the i^— integral now appears as an endpoint IR singularity as D ^ 4 
in the radial part of the i^^-integral. This is not surprising since we note from (|5.3|) 
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that the colhnear singularity of ( |5.2| ) is encountered on the plane spanned by the 
two massless momenta P and Q. The region of integration which contributes to this 
singularity is therefore identified by the limit K± — > 0. 

The combination of Sudakov methods and dimensional regularisation to study 



discontinuities in zero temperature type integrals is discussed in ||21|. In the present 
case though we have obtained a simplified structure for ( ^l2|) in terms of (|5l5| ), the 
dependence of the statistical factor on ko is still a complication. We therefore isolate 
the IR singularity in (|5.5|) into a simpler integral by rewriting 



Nko eiko) = [Nk, e{ko) - N^, e{xq)] + N^g e{xq). (5.7) 

The reason for this rearrangement is because now the term in square-brackets on the 
right-hand-side above gives an UV and IR finite contribution to ( p.5| ) ( and hence 
to Ji) as D ^ 4, while the simpler last term in (|5.7|) will give the pole ( and some 
finite parts). Calling the net contribution of the square bracket in ( ^.71) to Ji as Jia, 
we obtain (see Appendix B) 

- 1^- (5.8) 
The last term on the right-hand-side of (|5.7|) contributes to (|5.5|) the amount 



oo 



oo 



' '(-^i)^ + (ax + 7)' 



Since is a space-like vector with D — 2 independent components we may write 



K± = J2i=i Cj, where the Cj are orthonormal space-like vectors : Cj ■ ej = —6: 
Then (|0| ) becomes 



-00 ^ i-oo (af + + al_2) (E crff^sx) + (ax + 7) 



uiD-l)! dxl drr^-' /^'^^^'^^^ (5.II) 
^-00 Jo ^ + {ax + 7) 



, dx dyy""^-'- -i^ ^ (5.12) 

2 Jo Jo ^ ^ 2^ + (^a; + 7) ^ ' 

^ '-i2s)^'^-^ dx / dz -. 5.13 

2 ^ ' Jo e'^i + l Jo z + {ax + -i) ^ ' 



*The reader who does not find this too apparent may find a "posteriori satisfaction by examining 
the simphfied form of the relevant equation in Appendix B. 
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In going from ( p.lO| ) to |5.11| ) we switched from the Cartesian coordinates ai to spher- 



ical coordinates and the function u{D — 1) as defined by (|A.25| ) is the result of doing 
the angular integrals while / dr is the remaining radial integral. When (|5.13|) is 
substitued back into Ji, we will need the sum 

E (-^) 7^ V = E (-^) ( 7 7T + 7 7t1 • (5-14) 

The change of variables z = y{x + 1) for the first term in ( |5.14| ) and z = y\x — l|for 
the second term decouples the x — z integral in (|5.13 ) into a product of an x-integral 



and a simple ?/-integral. The resulting ^/-integrals can be evaluated explicitly using 
eqns. ( [Al26| - [A.28| ) . The x- integrals are still too complicated to be evaluated in closed 



form but the their integrands may be expanded about D = 4 directly or after some 



integration by parts. In this way we obtain the net contribution, Jib, of (|5.13| ) to Ji 

as 



J,^ = l^dqq rig |^^^[Boxl] + go (2[Box2] + 6 ln(g) [Boxl]) + Ag,[Boxl]j+0{D-A) , 

(5.15) 



where 



poo ]\f /.I jY _ jY roo 

[Boxl] = - / dx^^ + / dx-p^ + / dxqN^g{N^g-l)ln{x-l) , (5.16) 



(1-x) 



[Box2] = -—Ng - / dx^^lnxix + 1) + / dx ' Inxil-x) 
6 JO X + 1 Jo (1 — a^) 

+ fdx(^^\nx{x-l) + giV.,(iV,., - 1) , 

(5.17) 



and 



1 

^0 = 7, 5.18 

^" 5127r4 ' ^ ' 

-27r2 + 7r2ln2-37r2ln7r + 18C'(2) 
= 25fc5 ■ (5.19) 

Evaluating the integrals above numerically we obtain 
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with 



poo 

Ti = dq qnq[Boxl] 
Jo 



(5.21) 



-0.7167667897 (5.22) 



Thus finally 

Jl = JiA + Jib 



1 /3.33e _ 0.7167667897\ ^ ^^^^^^ 



Before going on to J2, let us recapitulate the story of Ji. Our objective was to 
obtain the order 1/(D — 4) and 0{D — A)^ pieces of Ji as D ^ 4. We first used a tech- 
nique popular and effective at zero temperature, the Sudakov decomposition, which 
exposed the singularity structure of Ji in a more manageable form. However, unlike 
T = cases, further progress was hampered by the presence of the statistical factor 
which in the Sudakov basis obtains a complicated energy dependence. Our next step 
was to isolate the singularity of Ji into a simpler integral by adding and subtracting 
terms in the original integral (see (|5.7| )). The simpler part of Ji, which we called 
Jib, had a statistical factor which did not depend on the transversal momentum K± 
of the Sudakov decomposition and so we were able in this case to proceed further in 
extracting the singular and finite pieces of Jib a.s D —>■ 4. The other part of Ji, which 
was created by the shift (|5.7] ) and which we called Jia, was a complicated but finite 



object at -D = 4 so we could evaluate it numerically. 



•J2 

In the evaluation of H2 in Sect. 4 we required. 



J2 = J21 + J22, (5.24) 



where 



■^^^ = -/ E (-^)'5i(^,7), (5.25) 



(27r)3(^-i) 



(T7 = ltl 
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and 

J22= ——^^N,N,N,5^{K')5^m5^{P') ^ (^7)52(^,7), (5.26) 
and where we have defined 

The evaluation of J2 1 is similar to that of Ji discussed earlier. We find 

= (5.28) 
Notice however that unlike Ji there is no singular 1/{D — 4) contribution from J21. 
Consider next the i^-integral in J22 (|5.26|) : 

In the Sudakov basis this becomes 

^-(^ i^il + (.x + 7)' 

with fco again given by ( |5.6| ). As D — > 4, the IR singularity in this integral is 
now more severe than the case of ( ^.5|) so that the rearrangement (|5.71 ) by itself is not 
sufficient to simplify the integral. However what is eventually required in (|5.26| ) is 

(a7)[eqn.(5.30)], (5.31) 

cr7=±l 

and we see that the IR behaviour of ( |5.31| ) as — is similar to that of ( |5.5| ). 
Hence the split ( p?7|) can be used as before to obtain the pole and finite pieces for 
J 22- Since the analysis is very much the same as before (see however some comments 
at the end of Appendix B), we simply state the result : 

1 / 0.6420474257\ 
^- = ^(l-«0-^U3irj- (5.32) 

Thus 

J2 — J21 + J22 

1 ( 0.6420474257\ 
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6 Results to Fourth Order 



The pressure up to and including order then follows from Sect. 2 and eqns.(|3T 



— 1 + -N\ - —— + 



T4 45 V 4 / 288 Un \ 3 



'^^.(0.4056) - e'N^ f ^ + ° In I) , (6.1) 



7r6 ' ' \^ 7r6 67r2 x 288 ^ 

The coupling above is an implicit function of the renormalisation scale fi. One may 
choose /i = T so as to eliminate the logarithm at this and higher orders. Then the 
pressure is a function of e(T). In principle the value of e(T) may be determined by 
comparing the perturbative calculation of some other observable at super-high tem- 
perature T (where the electron mass is negligible) with its experimentally measured 
value. Alternatively one can use the renormalisation group to relate e(T) to the cou- 



pling at some other scale A. Perhaps the most instructive thing to do is to write (|6Tl 
in terms of the perturbative renormalization-group-invariant coupling, at the energy 
scale T, given by 

e^(T) = e^(l + 01n^)+O(e^). (6.2) 
Defining a (T) = (T) / Arr we arrive at 

P_ ^ ^(i + In) _ ^ ^(T)N ( a{T)N \"' 
T4 45 ^ 4 ^ 72 TT 9v^ V / 



+ 



0.658 ±0.006 
N 



- 0.757 ± 0.004^ ( ^^^^^ ^ + O (a{T) 



5/2 



(6.3) 



Notice the disappearance of the logarithm. We have also indicated in ( |6.3| ) esti- 
mates of numerical uncertainties due to the evaluation of some integrals by quadra- 
tures. 

7 Fifth Order 

The next correction to the pressure is of order and comes about by dressing the 
photon lines of the 3-loop diagrams. Its calculation is completely analagous to that of 
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the term reviewed in Sect. 2 and the reader is encouraged to re-read the discussion 
there and in the references quoted. Since it has already been discussed at length in 
0, here we only sketch another derivation using the "ring-summation" formula. 

Consider the static, renormalised, one-loop photon polarisation tensor Il^^{qo = 
0,q). From gauge invariance, Q^U'^'^^Q) = 0, one obtains njo(0,g) = while exphcit 
calculations yield %(0,g 0) = 0{e^q^) and noo(0,g ^ 0) = + 0{e^q^). Thus 
for diagrams Gi and G2 (Fig. 2), one deduces from the usual power counting that it 
is only necessary to dress one of the photon lines with the static one-loop electric 
polarisation tensor to get the contribution. The resulting dressed diagarms are of 
the form of Fig. 5 and are contained in the full ring sum which is summarised by the 
formula (see the third reference of ||11||), 



Pring = -\ j [dQ] Tr{ln(l - D{Q) n(Q)) + D{Q) n(Q)} , (7.1) 

where the trace, Tr, is over Lorentz indices, fl is the full self-energy and D{Q) the 
bare propagator. Define F{q) = noo(0,g). Then the restriction of ( fZAl) to the static 
electric sector (as mentioned earlier this is the only sector which will give the 
contribution) gives 



2 J (2vr)^-i t2 ^ \ <1 

Now write F{q) = F^{q) + F^(g) + where the superscripts refer to the loop order. 
Then diagrams like those of Fig. 5 are obtained by truncating ( |7.2| ) to the appropriate 
sector, 

fen r7 911 ^ ^ f V (-^^ n F^jq) (F^g))-^ 

[eq.(/.2)J ^ -J j^^^ — ^ [7.6) 

T r d^-'q F\q) F\q) 

2 J (27r)«-ig2 (^2 + ^1^^)) • y' ■ ) 

By scaling q = mx one deduces that for the contribution it is sufficient to set g = 
in the of eq.( [7.4|) . Thus the order e^N^^"^ contribution to the pressure is given by 



[eq.(7.4)] ^ — F^(0) = ^ — V?U -0 (7-5) 
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where F^{0) = iv? = and P2 denotes the pressure at two- loop order at 

chemical potential /ie [0. In ( [7l5|) we used the relation noo(0,0) = e^d'^P/dnl [|lT|. 



The result ( [r.6|) may also be obtained by direct calculation of the left-hand-side of 

At order we still have the e^N^^'^ contribution obtained by dressing G3 (Fig. 4a). 
Now the ring summation formula ( [7.2|) has to be truncated in the sector where only 
iterations of F^{q) occur, but for one of them the subleading momentum dependence 
is taken while the rest are at zero momentum : 



Stt 9g2 '^^ 



(7i 



_ -e^T^N^/^ 7 - 1 + ln(4/7r) e^T^N^''^ ln(T//i) 

87rV27 12^^^ ^ Sttv^ Ovr^ ' ^ 

As in the last section, the logarithm eventually disappears when the pressure is written 
in terms of the temperature-dependent coupling ( |6.2| ). 

The derivation of (|7.5| ) may be extended naturally to re-obtain the identity of 
linking the contribution of the pressure at order 2n + ?> {n > 1) from diagrams 
with one-fermion loop to the pressure at order 2n. The relation ( |7.5| ) and its higher 
order extensions were stated in p| for the case of massless fermions at zero chemical 
potential. Clearly one can relax these restrictions since only the dressing of the photon 
is involved. Thus in general 



e^TJF^iO) d^P^^ 
^"+3 " Stt ^^-'^^ 



relates the gauge- invariant pressure at order 2n + 3 (n > 1), from diagrams with one 
fermion loop, of QED with massive electrons at non-zero temperature but arbitrary 
chemical potential, to the pressure at order 2n (the nonzero T is required so as to 
isolate the zero mode to be dressed). In this general case, F^{0) is the lowest order 
(e^) electric screening mass at nonzero T, /^e and electron mass. 
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8 Conclusion 



With a(T) = e2(T)/47r, and defining g'^ = a{T)N/Tc, the pressure of QED with N 
massless Dirac fermions at nonzero temperature, T, is given to fifth order by 



P 

2^4 



ao + /as + g^as + g^a^ + h/N) + ^^(05 + h/N) + 0(/) , (8.1^ 



with 



57r2 



a2 



as 



72 

27r2 



9v^' 

04 = -0.757 ±0.004, 
64 = 0.658 ±0.006, 

9^3 



-7r2 



2^3 



-2.849... 



8.2) 
8.3) 

8.4) 

8.5) 
8.6) 

8.7) 
8.8) 



Real world QED corresponds to = 1, but since we have ignored the electron 
mass the results are applicable only at extremely high temperatures Numerically, 
the fourth and fifth order terms we have found are small corrections in the regime 
where the coupling itself is small. However, since perturbative QED is not asymp- 
totically free, the effective coupling a{T) increases slowly with temperature so the 
results might be of use for physics of the very early universe or, more speculatively, 
for studies of strongly coupled QED (some references are in 0]). 

It is an amusing fact that the order contribution to the pressure of QED is 
much easier to calculate than the contribution. Indeed, as discussed in the last 
section, the fifth order calculation may even be extended to massive electrons and 
nonzero chemical potential (but nonzero T). However we have not bothered to give 
the explicit expressions in those cases because the fourth order calculation at nonzero 
T is itself only known for massless electrons at zero chemical potential 

A 3-loop calculation in QCD will differ from the QED case in two respects. Firstly 
there is an increase in the number of diagrams. This however is not a problem (except 
for tedium) as we feel that our approach using the frequency-sum algorithm discussed 
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in Sect. 5 and the Sudakov method of Sect. 6 is general enough to handle any new in- 
tegrals that might arise. The second difference is that the static electric polarisation 
tensor in QCD behaves as noo(0, 5' — > 0) = + qT and this gives rise to the g'^ \ng 
term ^ from the sum of ring diagrams. Thus in this case one has to be more careful 
in using dimensional regularisation (as in this paper) to extract this term and also 
the constant under the logarithm. 

Note added in proof: The three-loop free energy of hot Yang-Mills theory has been 
obtained in Ref.[^ . 
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Appendix A 

Al. Identities for statistical factors 

Let Tlx = l/{e^ — 1) and = l/{e^ + 1). Also, denote the step function by 9{x) and 
the sign function by e(x). We have 



n-x 


= -{l + Ux), 


(A.i; 


N-x 


= l-Nx, 


(A.2; 




= -9{-x) + e{x)n\x\, 


(A.3; 


Nx 


= 9{-x) + e{x)N\x\, 


(A.4; 






(A.s; 


nx+y{N_x - Ny) 


= NxNy, 


(A.e; 


Nx+y{N^x + riy) 


= N n 


(A.7; 


^x+y+z{^x^z + n^y^-z ~^ ^x^^-y) 


TT'x^y^Z) 


(A.s; 


Nx+y+z{N_xN_y - N_xN, + NyN,) 


= NxNyN,. 


(A.9; 



The equations (|A.1|) - (|A.7|) follow from the definitions of Ux and Nx while the last two 
are obtained by iterating ( [A.5D -( [A77D . 



A2. Standard Results 



00 -1 

E — 



n=l 



E 

n=l 



71+1 



Mi{a) 
M2ioi) 



dx- 



X 



— 1 

dx 

+ 1 

00 x^~^ 
dx 



C(«), 

1 - 2^-") C(«) , 
r(a + l)C(a + l), 
(l-2-)r(« + l)C(a + l) 



1 „ _ ^ 



(AlO) 
(All) 
(A12) 
(A13) 



— y 



—Tip) V^, 

{p> 0, r > 0, -1 < y < 1), (A.14) 
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rco 1 

/ dxx^'-^e-^'' = —T(u), ReU.v) > 0. (A.15) 

Jo fJ," 

^dxx^-\l-xy-^{ax + b{l-x)+cy^^^''^ = (a + c)-'^(6 + c)-"5(/i,z/), 

a>0, 6>0c>0, 



Re ii> 0, Re u > 0. 

(A.16) 

A3. Derived Relations 

/•oo roo 

M3{a,(3)= / dx dy x" n^+y = r{a + 1) r{(3 + 1) ({a + (3 + 2) , (A.17) 
JO Jo 

Al4(a, P)= dxdy x" y^N^+y = (l - 2-(^+"+^)) r(a + 1) r(/? + 1) C(a + + 2) 

(A.18) 

The result (|A.17|) is obtained by using in sequence the relations ( |A.14|) , ( |A.15| ) and 
( lAlOl ), and similarly for ([Alsl) . 



• Some Covariant Integrals. 
Define 



= ^ K.{P + Q)±P.Q {K.Qr^ ^^-^^^ 

= ^ K.{P-Q)±P.Q{K.Q)^- ^^-^^^ 



Then 



F^(+,±) = 5+(Q2)5^(p2W_J^ C^(+,±), (A.21) 



F^(-,±) = 5^(g2)5^(p2) C^(-,±), (A.22) 



where 



uj{D - 1) 2^-5-w 

giv(-,±) 
cj(D - 1) 2^-5-^ 



5(D/2 - 1, D/2 - 1 - A^) Pn{±), (A.23) 

5(3 + A^ - D, D/2-1- N) Pn{±) (A.24) 
_5(3 + iV-D, D/2-1) P^iT), 
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and 



uj{D) 

PNi + ) 



D-l 

2n— 



r(s=i; 

TT csc((A^ + 3 -D)7r) 



(A.25) 
(A.26) 



TT cot {{N + 3- D)-k) 



(A.27) 



The result ( |A.21|) is obtained as follows : Since the integral is covariant, it may 
be evaluated in any convenient frame. Choose p = —q. The the only nontrivial 
integarls are the radial integral dk and the angular integral j\d cos 6 , where 9 
is the angle between k and q . These two integrals may be decoupled by a simple 
change of variables and the angular integral then evaluated using ( |A.16D while the 
radial integrals are 



oo ^D-3-N 

dz — — — 

z ±1 



(A.28) 



For the + case the integral is standard 27 1 while the — case is interpreted in the 



principal value sense and the result is indicated in (|A.27|) . Consider instead the 
integral 



PNi-) = 



dz 

-2 - 1 =F ^0+ 



(A.29) 



(A.30) 



Compared to P]\}{~), the integral Pn{—) has the original pole at z = 1 shifted above 
or below the real-axis and this shift may be viewed as a "regularisation" of P/v(— ) 
which is then given by the real part of Pn{—)- Combining eqns. ( |A.26| - |AI27| ) and 
( |A.3CI| ), we can relate -P/v(— ) to Pn{+) which is regular , 



N{ 



(A.3i: 



A direct derivation of this result can be obtained by an analytic continuation of P/v(±) 
for 2; < 0. We omit this presentation since it requires more involved considerations 
of the integrals for 2; < 0. 
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A4. Relations between Gamma and Zeta functions 

One has the standard formulae [^| 



f-ir /111 \ 

n-n + e) = - + (1+ +...-- 7) + 0(6)) (A.32) 
7 = lim (c{z) —) = 0.5772157... , (A.33) 



2^1 V z-1 
7i'-X{z) = r r(l - z) C(l - z) sin — , (A.34) 



and one may deduce 

C(0) = C'(0) = -^ln 27r, C"(0) = -2.00635645..., 

C(-l) = C'(-l) = -0.16542114369..., 

C"(-l) = -0.2502044...., 

((2) = — , C'(2) = -0.937548254..., 
6 

C"(2) = 1.9892802342... 

The specific values above may be obtained by a Taylor expansion of both sides of 
( |A.12| ) and ( |A.13| ) with respect to an appropriate value of a, and using ( |A.32| )-( |A.34| ). 



Appendix B 

Here we describe how to obtain ( |5.8| ). The contribution of the square bracket in 
to ( |5.5|) is finite as D 4 and is given by 

v^Ti- l<eK., ' . I - ^^'^ (B.i) 



oo 



2xs 



where s = P- Q=pq — p-q and ko = -^^p + xg + {K^_)q. As in the discussion 
of ( p.9|) we first decompose along two orthonormal space-like vectors. Thus we 
write 

= cTie^i) + ^26^2) (B.2) 

where the are two orthonormal space-like vectors : e^^-* ■ e^^-* = and (e*^^-')^ = 
(^g(2)^2 _ However unlike the case of (|5.9|) the integral in ([B.l|) depends explicitly 



on the complicated energy /cq and so we now require an explicit parametrisation of 
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our basis in order to proceed. Though the Lorentz symmtery of the integrals in Ji, of 
which ( |B.1|) is a part of, is broken by the heat bath, we still have three dimensional 
rotational invariance in the {p, q) integrals. Therefore we can choose the following 
explicit basis (we are grateful to Cosmas Zachos for discussions on this point) for the 
evaluation of ( |B.1| ) and its contribution to Ji : 

P^ = p(l, 0,0,1) 

= g(l,0,sin 0,cos 0) 
e(i) = (0,l,0,0) 

e(^) = (T^,0,l,-^). (B.3) 

1 — cos (p 1 — COS (p 

where is the angle between p and q. 

The basis ( p.3| ) satisfies the requirements 

P^ = Q^ = e(i) ■ P = e(i) ■ Q = e(2) ■ P = e^^) . Q = , (B.4) 



and Po = P, qo = q,s = P- Q= pq{l — cos 0) > (— vr < < vr). Thus ( |B.1|) can 
be written as 



oo roo roo 



1 + ^2 ^ ^' + {ax + 7) 



00 J —00 J —00 

2xs 



with fco = -^^^^P + xq + a2e[)^\ 

Now it is convenient to transform from the Cartesian {(Ti,a2) coordinates to the 
polar coordinates defined by 

ai = y sin 9 , (^2 = y cos 9 

and then change variables y"^ ^ y to obtain 

v= r dx r^rd9 1^^- ^(^q)-^-^ ^(^g)] me) 

J-00 Jo y Jo 2^ + (^a; + 7) 

with k'Q = ^ + xq + y/ycos 9 ey . 

Next split the x— integral in ( [B.6| ) as /q°° dx + dx and in the second part do 
X — >• —X, 9 ^ IT — 9, and note that eventually we need J2'y=±i V, to get (after some 
simplification) 



7=±1 7=±1 



°° . dz 



Y.V=Y: / dx - d9 ^ (B.7) 



z Jo z + {ax + 7) 
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where 



pz + xq + V2sxz cos 6 e^'^ 

= pz + xq + J2pqxz{l + cos 0) cos 9. 



(B.8) 



Furthermore it follows from ( p.8|) that ko > 0, and thus we can write (|B.7|) as 



7=±1 7=±1 



, /"^ dz 
dx — 



d9 



z Jo z + {ax + 7) 



(B.9) 



Finally Jia is given by 

JlA = 



d^p d^q {Np + Up) Ug 



Apq 

(47r)(27r) 
4(27r)9 



7=±1 



(27r)9 
dp p {Np + Up 



POO pi 

I dq q Uq I d cos (j) V , (B.IO) 

^0 j-i til 



with V defined through equations ( |B.9| ) and ( |B.8D . The integrals in ( [B.1CI| ) can now 
be performed numerically. One technical point that should be noted is the principal 



value singularity in the factor 



1 



occurring in ( p.9|) and ( p.lOp . Consider, for 



example, the case Pjz^zzt) ( l^-^D - Doing the change of variables z = y\x — l\ gives 
us the factor to deal with in the new /q°° dy integral. We use the definition of 

the principal value to write 



1 



lim 

<5^0 



dy P — ^ X (rest) 
I-'' dy 



dy , . dy , , 

/ -(rest) + / -(rest) 

JO y — i- y — i- 



:b.ii] 



where "(rest)" denotes the rest of the integrand in ( p.9|) after the change z = y\x — l\. 
In the second term in (|B.11|) one can do the change of variables y — > ^ and then 
combine the result with the first term to get an integral of the form 

(new rest) , (B.12) 

Jo y-l 

where the limit 6^0 can be taken because now the integral is finite as y 1 since 
(new rest) — > (y — 1) as y — * 1. Once the principal value singularity has been removed 
by transforming ( [B.ll| ) into ( |B.12| ), the numerical integration of ( [B.lOp can proceed 
and we obtain 

= (B.13) 



2%^ 
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Some remarks are in order. The change of variables z = y\x — 1| that yielded 
( p.ll|) from ( p.lO|) results in an x-integral of the form /q°° dx{-)/\x — 1|. In this case 
one appears to have created a new singularity at x = 1. However this is not true 
because the rest of the integrand actually compensates it as s ^ 1. In other places 
in this paper, for example in the evaluation of Jib ( |5.14| ) and also -F/v(— , ±) ( |A.20| ) 
near D = 4 dimensions, the change z = y\x — 1\ actually creates a singularity at 
X = 1 when = 4. Noting that dx{-)/\x — 1| actually means Jq dx{-)/{l — x) + 

dx{-) / {x — 1), we have used dimensional continuation separately for each part to 
regulate the singularity at x = 1 in those cases. 

For another example, consider an integral of the form 



^ / dx 



1 



z — (x — V 



+ 



1 



z+ {x- 1) 



(B.14) 



which appears in the evaluation of J2 ( |5.31| ). If one now naively does the change 



z = y\x — 1| to simplify the integrals, the x-integral becomes /q°° j^^j^^ (new*) and the 
singularity at x = 1 is not compensated. The safe way to proceed in this case is to 
first do X ^ xz in ( |B.14| ) to get 



dx 



°° dz 



1 



z(l-x) + l) 



+ 



1 



z{x + 1) 



X 



xz) . 



(B.15) 



Then do 2; — * z\x — 1| for the first term, z z{x + 1) for the second, and finally 
proceed as per ( |B.11| )-( |B.12| ) to remove the principal value prescription so that the 
integrals may be handled numerically. 
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Figure Captions 

Fig.l : 

Contribution to order e^. The wavy line represents the photon propagator. 
Fig.2 : 

The order e'^N contributions Gi and G2- 
Fig.3 : 

Ultraviolet counterterm diagrams for Fig.2. 
Fig.4 : 

Fig. 4a is the e^A^^ contribution (G3) while Fig.4b is the corresponding counterterm 
diagram Xi. 
Fig.5 : 

Diagrams obtained by self-energy insertions along one photon propagator in the dia- 
grams of Fig.2. 
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This figure "figl-l.png" is available in "png" format from: 
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This figure "fig2-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9409269vl 



